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Abstract 

We report measurements of turbulent drag reduction of two different 
polyelectrolyte solutions: DNA and hydrolyzed Poly aery lamide. Chang- 
ing the salt concentration in the solutions allows us to change the flexibility 
of the polymer chains. For both polymers the amount of drag reduction 
was found to increase with the flexibility. Rheological studies reveal that 
the elongational viscosity of the solutions increases simultaneously. Hence 
we conclude that the elongational viscosity is the pertinent macroscopic 
quantity to describe the ability of a polymer to cause turbulent drag re- 
duction. 

The phenomenon of turbulent drag reduction describes the diminution of the 
dissipation in turbulent flow by adding tiny amounts of polymers. In this way for 
instance the amount of liquid that is transported in a pipe for a given pressure 
drop can be increased significantly. This implies a wide range of industrial 
applications, and consequently the effect, discovered in 1949 by Toms P, has 
been studied intensively over the past decades (see e.g. Ref. |2]). In spite of this 
effort, the underlying physical mechanisms are ill understood and none of the 
existing theories match the existing experimental data [3] . The key macroscopic 
property of these solutions that is known to be significantly different from that 
of the solvent is its elongational viscosity. Therefore, drag reduction is usually 
attributed to the elongational viscosity; however the precise connection remains 
unclear. 

The elongational viscosity describes the resistance of a liquid to an elonga- 
tional flow. In Newtonian liquids the elongational viscosity is just given by the 
shear viscosity times a geometrical factor of three, the Trouton ratio W. The 
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addition of a small amount of polymers to the solvent does not change the shear 
viscosity significantly, however the elongational viscosity can increase by several 
orders of magnitude. Taking an elevated elongational viscosity into account, a 
qualitative explanation of drag reduction has been given by Landahl [S]: He sug- 
gests the elongational viscosity can suppress the occurrence of streaks - regions 
of flow with a high elongation rate - in the turbulent boundary layer. These 
boundary layer instabilities lead to "blobs" of fluid that are ejected from the 
boundary layer into the bulk of the liquid and which generate the turbulence. 
An observation that provides support for the idea that the boundary layer is 
important was made by Cadot et al. 0. They showed experimentally that 
drag reduction occurs only in boundary layer driven turbulence: if no bound- 
ary layers are present, indeed no drag reduction is observed upon addition of 
polymers. This implies that it is the way the turbulence is generated by the 
boundary layers that is altered by the addition of polymers. 

To our knowledge, however, an explicit experimental proof of the relation 
between the elongational viscosity of a given polymer solution and its ability to 
cause turbulent drag reduction is still missing mainly because the determination 
of the elongational viscosity of drag-reducing polymer solutions has proven diffi- 
cult. It is the aim of this publication to provide this relation. In order to do so, 
we study both the elongational viscosity and turbulent drag reduction of poly- 
mers that allow us to tune their chain flexibility. The large elongational viscosity 
of dilute polymer solutions is usually attributed to the resistance to stretching 
of the polymer chains in the elongational flow field. Therefore, flexible poly- 
mers should have a higher elongational viscosity. If in addition the phenomenon 
of drag reduction is indeed related to the elevated elongational viscosity, drag 
reduction should also increase with increasing chain flexibility. We show here 
that this is indeed the case using polyelectrolytes, for which the chain flexibility 
can be tuned by the addition of salt to the solvent [?]• In that way the chain 
flexibility is altered without changing either the polymer chemistry or the chain 
length (the effects of which on drag reduction are still ill understood) . 

The first polyelectrolyte solution we use are aqueous solutions of double 
stranded A DNA (48.5 bp) at the overlap concentration c* of AO^g/ml. DNA 
is a model monodisperse polymer and its physical properties, for instance the 
relaxation time, are very well defined 0. Our solutions contained a lOmM Tris 
and ImM EDTA buffer (alternatively we also use a 2mM Tris buffer) and the 
NaCl concentration was varied from to lOmM, changing the persistence length 
by a factor of three and the relaxation time from about 10 to 30ms. The second 
polyelectrolyte used was AQ^g/ml hydrolyzed Polyacrylyamide (HPAA) with a 
molecular weight of 5 x 10^ amu in water with NaCl concentrations varying 
from ImM to 18mM. 

First the shear rate dependent viscosity rj (7) of our samples was measured by 
use of a standard rheometer (Reologica Stress Tech) with a cone plate geometry 
that assures a laminar fiow (figure 1). In the range 20s~^ < 7 < 2000s~"'^ 
all HPAA samples show shear thinning, the effect becoming less pronounced 
at higher salt concentrations. This was to be expected, as stiffer polymers 
show in general stronger shear thinning. The DNA solutions showed only a 
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slight shear thinning and no dependence on the salt concentration could be 
detected. At at a shear rate of 7 = 2000s~^, comparable with the highest 
shear rates used in the drag reduction experiments. The shear viscosities for 
the aqueous HPAA solutions varies in between 1.22mPas < 77 < 1.38mPas 
with the salt concentration and for the aqueous DNA solutions the viscosity is 
77 ~ l.OSmPas. Our rheological measurements did not reveal any measurable 
normal stresses within the resolution of iVi ~ IPa of the rheometer for the 
aqueous AO^g/ml solutions employed here. 

To investigate turbulent drag reduction we measured the drag of a turbulent 
flow of the different liquids in a Couette Cell with a gap oi 6 — 1mm. Even 
if the industrial implications of turbulent drag reduction are more likely to be 
pipe flow, it is experimentally much more convenient to study Couette flow; 
this is in fact a rather common way to characterize polymer solutions for their 
ability to exhibit drag reduction A further advantage of this setup with a 
high surface to volume ratio is to probe directly boundary layer effects. 

Drag reduction measurements were performed on our standard rheometer. 
It was driven in controlled stress mode, and drag reduction DR is defined as 
the normalized difference of the drag d on the inner cylinder by pure water 
and by the polymer solution respectively at same Reynolds number DR = 
{dH20 — dpoi) /rfffaO X 100%. The Reynolds number is defined as Re = 
being the rotation speed of the inner cylinder, p the liquid density and -q the 
solvent viscosity. 

Figure 2 shows the measured drag reduction as a function of time for a 
Reynolds number Re w 2000. Both for the DNA and for the HPAA solution 
DR is found to depend on the salt concentration in the solvent: drag reduction 
increases with increasing salt. For the DNA solution the results are qualitatively 
as for the HPAA, but the effects are less pronounced. 

We believe the interesting temporal behavior for times t ^< 400 sec to be a 
rather complicated mixture of disentanglement of different polymer chains [S], 
rupture of the individual polymer chains by the strong flow [5] , and relaxation 
of the flow towards a steady state. In addition to the data shown here, at longer 
times we observe a tendency to smaller values of drag reduction, probably due 
to further degradation of the polymers [H] . 

Drag reduction as a function of the Reynolds number is shown in figure 3. 
For each new Reynolds number the system was given 60 seconds to reach a 
steady state, a reasonable compromise in view of the different temporal effects 
discussed above. For all the solutions, we observe an increase of Drag reduction 
with increasing ionic strength. We checked that this is not due to our use of 
the Couette cell, or the very moderate Reynolds numbers employed here by 
repeating the experiments in the turbulence cell of Cadot et al. [HI- These 
high-Reynolds number experiments gave very similar results. The highest Drag 
reduction measured with HPAA for a Re « 4.10^ was DR ~ 30% rather close 
to the drag reduction DR ^ 20% measured with the Couette cell. The complete 
study of drag reduction at high Reynolds number will be published elsewhere 

cni 
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Having measured the effect of salt on drag reduction^, we now turn to the 
elongational viscosity. As the velocity fluctuations in turbulence are violent, 
and rapidly change direction, it is unlikely that the polymer extension in such 
flows reaches a steady state. The pertinent experiment is therefore to study 
incipient (startup) elongational flow, and measure the response of the polymer 
solutions to that . 

In order to measure the dynamic elongational viscosity 'i]c,{t) of our solu- 
tions we followed the method proposed by Bazilevskii et al. |T2] and further 
developed by Amarouchene et al. Therein the authors describe the de- 

tachment process of a droplet of a diluted polymer solution from a capillary. 
The addition of a small amounts of flexible polymers to water inhibits the fi- 
nite time singularity break up ^31 of the droplet, and a cylindrical filament is 
formed between the capillary and the droplet. The flow profile in the filament is 
purely elongational and by balancing the capillary forces to the elastic stresses, 
the elongational viscosity 77e(i) can be extracted from the measured filament 
diameter h{t). The elongational viscosity i]e{t) is time dependent because in 
the course of the experiment the polymers are more and more stretched by the 
flow. The experiments reveal that the temporal behavior of h{t) is exponential, 
implying a unique thinning regime for which the elongational rate e = 
remains constant. 

Recently Anna and Mc Kinley have shown that filament thinning tech- 
niques analogous to the one used here lead to measures of the elongational viscos- 
ity that are consistent with the more established filament stretching techniques. 
The latter cannot be employed here since it only allows for measurements for 
much more viscous solvents (typically IPa.s) |15j . 

It turns out that for the aqueous solutions used in the DR experiments, the 
droplet detachment process is too fast to observe, even using a rapid camera 
(Photonetics) working at a 1000 frames/s. Increasing the solvent viscosity slows 
down the dynamics of the filaments sufficiently so that these can be followed, 
and the elongational viscosity can be determined. We will show below that these 
measurements suffice also to deduce the elongational viscosity of the aqueous 
solutions used in the drag reduction experiments. 

The solvent viscosity could be adjusted by adding a certain amount of glyc- 
erol. Figure 4 shows the photographs of the detachment process for the DNA 
solution for the OmM and the lOmAI salt concentration. The trend can be 
observed directly on the pictures: the more flexible polymer has a higher elon- 
gational viscosity; it is evident from the figure that the filament thinning is 
much slower. 

In the filaments, for the DNA solutions, the elongational rate varies with 
NaCl : e = 90s~^ for low concentration and e = 50s^^ for the high concentra- 
tion. In order to be able to compare different experiments, it is thus useful to 
introduce the so called Hencky strain, the product of the elongation rate and 

^In a recent paper, Choi et al. 1111 have reported on the drag reduction properties of 
DNA at lower concentrations but higher Reynolds numbers . However, their focus was on 
studying the degradation introduced by the high stretching fields and no correlation with the 
elongational viscosity was given. 
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the elapsed time ei, which measures the total deformation of the fluid element 
containing the polymer. We define the zero of the Hencky strain for the moment 
when the filament is formed, i.e. when the flow is strong enough to interact with 
the polymers 

Figure 5 shows the elongational viscosities of the different samples, for dif- 
ferent salinities as a function of the Hencky strain. Again both for the DNA 
and the HPAA, rje depends strongly on the salt. For HPAA, rje changes more 
than a factor of three between the lowest and the highest salinity, while there 
is still a factor of almost two observed for the DNA solution. 

A prediction of all polymer flow models is that the elongational viscosity 
scales linearly with the solvent viscosity ^Hl- If this is true, we can directly 
extract the elongational viscosity of the aqueous solutions used in the drag 
reduction experiments from the elongational viscosities of the glycerol- water so- 
lutions. To verify this, we prepared a series of additional samples of 300fig/ml 
PAA in different glycerol-water mixtures. This allows us to change the solvent 
viscosity 77 in the range lOmPas < SOmPas. The inset of figure 5 shows the 
Trouton ratio Tr, the measured elongational viscosity r]e of the polymer solu- 
tions divided by three times the viscosity of the solvent (Tr/'i = 1 for Newtonian 
liquids). The Trouton ratio of the different solutions perfectly collapses on to 
a single master curve. For even smaller viscosities slight discrepancies can be 
observed, but here the dynamics becomes again too fast to allow for a precise 
determination of rje- The perfect collapse also indicates that the chain flexibility 
is not modified by the addition of different amounts of glycerol. Therefore we 
can obtain the elongational viscosities rie of the aqueous polymer solutions by 
dividing the data from the glycerol water solutions by the ratio of the solvent 
shear viscosities. 

It is therefore now possible to relate the drag reduction directly to the elon- 
gational viscosity. The result is shown in figure 6. We plot the drag reduction 
for a given Reynolds number as a function of the elongational viscosity rj^ at 
a Hencky strain of 1. Drag reduction is shown to increase monotonically with 
the elongational viscosity for a given polymer and the data for the different 
polymers collapse, demonstrating that indeed the elongational viscosity is the 
pertinent macroscopic quantity to account for DR. As far as we can tell, this 
is the first explicit demonstration that the turbulent drag reduction is directly 
connected with the elongational viscosity of a polymer solution. 

These results, especially those on DNA, open the way to a microscopic un- 
derstanding of both the enormous elongational viscosity dilute polymer solutions 
can have, and the surprising results thereof: turbulent drag reduction. It has 
recently become possible to observe the deformation of the polymer chains due 
to flow of single DNA molecules using fluorescence microscopy |H|. Perform- 
ing such an experiment on the filaments would allow to relate the extension of 
the individual polymer chains, to the macroscopic stresses the chain extensions 
generate. If subsequently the polymer conformation is visualized in a turbulent 
boundary layer, an estimate of the extra stresses in the boundary layer could 
be obtained, which could be the key to the understanding of this phenomenon. 
These experiments are in progress. 
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figure 1: The shear dependent viscosity of the different polymer solutions. 
Filled symbols: aqueous HPAA solutions ImM NaCl (circles), 2mM(up trian- 
gles), 7mM (down triangles), 18mM(diamonds). Filled symbols: aqueous DNA 
solution OmM NaCl (up triangles). 

figure 2: The turbulent drag reduction DR as function of time at a constant 

Stress of 8 Pa yielding a Reynolds Number Re w 2000. Filled symbols: DNA 
solutions lOmM NaCl (down triangles), the other symbols like in figure 1. The 
error bars indicate the statistics of three consecutive measurements. 

figure 3: The turbulent drag reduction DR as a function of the Reynolds 
number Re for 4QjjLg/ml HPAA a) and AQug/ml DNA solution. 

figure 4: The detachment process of the droplet for the AO fig /ml DNA so- 
lutions with %Q/2Qvol% glycerol- water as solvent and a),b),c),d) OmM Nacl and 
e),f),g),h) IOtoM NaCl. The rows show different time steps t in regard to the 
moment tc when the filament breaks. a),e): tc — 120ms; b),f): tc — 94ms; c),g): 
tc — 53ms; d),h): tc — 20ms. The viscosity of the solvent is 80 mPas and the 
diameter of the capillary used to study the droplet formation has a diameter of 
1.2mm. 

figure 5: The elongational viscosity 77,. for the different polymer solutions. 
The symbols are like in figure 1, only the solvent is a 80/20vol% glycerol- water 
mixture. Inset: The Trouton ratio Tr for solutions with SOO/i^/mZ PAA in dif- 
ferent water glycerol mixtures and OmM salt. The solvent viscosities are 7]= 80 
(circles), 43 (squares), 20 (up triangles), 10 (diamonds) mPas. 

figure 6: The drag reduction DR at a Reynolds number Rcpoiy = 1400 as 
a function of the elongational viscosities of the aqueous polymer solutions, for 
different salinities at a Hencky strain oi et = 1. To allow for a comparison be- 
tween samples with a different shear viscosities, the Reynolds number Rep^iy is 
calculated using the laminar shear viscosities of the polymer solutions at a shear 
rate 7 = 2000. Filled squares: DNA-solutions, open squares: HPAA solutions. 
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